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Abstract 

We report on the first measurement of the single spin analyzing power (An) at -y/s = 200 GeV, 
obtained by the pp2pp experiment using polarized proton beams at the Relativistic Heavy Ion 
Collider (RHIC). Data points were measured in the four momentum transfer t range 0.01 < |t| < 
0.03 (GeV/c)2. Our result, averaged over the whole |t|-interval is about one standard deviation 
above the calculation, which uses interference between electromagnetic spin-flip amplitude and 
hadronic non-flip amplitude, the source of A^. The difference could be explained by an additional 
contribution of a hadronic spin-flip amplitude to A^. 

PACS numbers: 13.85.Dz and 13.88.+e 
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I. INTRODUCTION 



The pp2pp experiment |l[ O, S] at RHIC is designed to systematically study polarized 
proton-proton (pp) elastic scattering from ^/s = 60 GeV to ^/s = 500 GeV, covering the 
range from the region of Coulomb Nuclear Interference (CNI) to 1.5 (GeV/c)^. Studies of 
spin dependence of pp scattering at small momentum transfers and at the highest energies 
presently available at RHIC offer an opportunity to reveal important information on the 
nature of exchanged mediators of the interaction, the Pomeron and the hypothetical Odderon 
(see Ref. 0,0] and references therein). The theoretical treatment of small-t scattering is still 
being developed, hence the experimental data are expected to provide significant constraints 
for various theoretical approaches and models (see Ref. and references therein). 

In this paper we present the first measurement of the analyzing power A^- in pp elastic 
scattering of polarized protons at RHIC at ^ = 200 GeV and 0.01 < |t| < 0.03 (GeV/c)^ 
An is defined as the left-right cross section asymmetry with respect to the transversely 
polarized proton beam. In this range of t, An originates mainly from the interference 
between electromagnetic (Coulomb) spin-flip and hadronic (nuclear) nonflip amplitudes Qj. 
However, it was realized that in the Coulomb-nuclear interference (CNI) region is a 
sensitive probe of the hadronic spin-flip amplitude 0. A possible hadronic single spin-flip 
amplitude would alter An and its effect would depend on the ratio of the single spin-flip 
amplitude (^5) to nonflip amplitudes (0i and ^3), Eq.([H): 

rs = m05/(v^ lm(0i + 03)/2), (1) 

where m is the nucleon mass (see Ref. jiSii] for definitions). 

Other measurements of An performed at small t have been obtained at significantly 
lower energies, by at least a factor of 10, than the present experiment. These measurements 



include recent hi 
for elastic pp jS 



iprecision resu 



ts from the RHIC polarimeters obtained at i/i = 13.7 GeV 
' lo| scattering, as well as earlier results from BNL AGS 



c Land pC |8| 

for pC scattering h| at a/s = 6.4 GeV and from FNAL E704 for pp scattering |l^ at 
v/i = 19.7 GeV. 

The combined analysis of the present result with the earlier ones, especially with the 
very accurate results ot Rets Q, will help to disentangle contributions ot various ex- 
change mechanisms involved in elastic scattering in the forward region ^3|]. In particular, 



such analysis will allow us to extract information on the spin dependence of the diffractive 
mechanism dominating at high energies. 



II. THE EXPERIMENT 



The two protons collide at the interaction point (IP), and since the scattering angles 
are small, scattered protons stay within the beam pipe of the accelerator. They follow 
trajectories determined by the accelerator magnets until they reach the detectors, which 
measure the x, y coordinates in the plane perpendicular to the beam axis. Those coordinates 
are measured by Si detectors in the Roman Pots, which are positioned at the location 
that satisfy so called "parallel to point focusing". More details on the experiment and the 
technique used can be found in [3, B]- The layout of the experiment is shown in Fig. 
The identification of elastic events is based on the collinearity criterion, hence it requires the 
simultaneous detection of the scattered protons in the pair of Roman Pot (RP) detectors Q| 
on either side of the IP. 

The elastic event trigger required a coincidence between signals in the RP's scintillators, 
belonging either to arm A or arm B, see Fig. ^ For each arm the trigger counters in RPl 
and RP3 were used. The overall trigger was the logical OR of a coincidence between up (U) 
and down (D) pots: (RP3U AND RPID) OR (RP3D AND RPIU) in coincidence with the 
beam crossing signal derived from the RHIC master clock. 



III. SELECTION OF ELASTIC EVENTS 



The detectors in the inner Roman Pots were used for elastic event reconstruction, as this 
provided the highest acceptance for the experiment. Particle hits in the silicon detector 
were identified for each strip requiring that the energy deposited (AE') was /\E > 5a of its 
pedestal value. From those hits a cluster of consecutive strips was formed and the coordinate 
for that cluster was calculated as an energy-weighted average of the positions of the strips. 

For each RP a hit was formed for an (x,y) coordinate using the clusters in two x planes 
(xi, X2) and two y planes {yi, y2)- A hit required that the distance between two clusters 
from adjacent planes was |xi — X2I < 2 strips, the same for y-coordinate — 1/2I < 2 strips. 
For matched clusters a single x and y coordinate was calculated as an arithmetic average of 
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FIG. 1: Layout of the PP2PP experiment. Note the detector pairs RPl, RP2 and RP3, RP4 he 
in different RHIC rings. Scattering is detected in either one of two arms: Arm A is formed from 
RP3U and RPID. Conversely, Arm B is formed from RP3D and RPIU. The coordinate system is 
also shown. 

the two. In case there was no match with the second plane one coordinate was used. 

Because of the coUinearity of the scattered protons one has to require a correlation be- 
tween coordinates measured on each side of the IP. Hence the main criterion to select the 
elastic scattering events was the hit coordinate correlation in the corresponding silicon de- 
tectors on the opposite sides of the IP. An example of the correlation of the x-coordinates 
of the detected protons is shown in Fig. |21 Note the diagonal band of the elastic events and 
relatively small background. 

Since the events for which the protons were detected in all four RP's allowed reconstruct- 
ing of the momentum vectors of the scattered protons at the detection point, a subset of 
those events was used to get better knowledge of the mean coordinates of the collision vertex 
and of the mean angles of the beams in the IP. The mean values and widths of those distri- 
butions were also used to determine the correction to the calculated transport matrices, and 
the beam position at detectors in the horizontal plane. The widths of these distributions are 
dominated by the beam emittance of about IStt mm ■ mrad for both x- and y-coordinates 
and by an uncertainty of about 60 cm (rms) of the vertex position along the beam axis. 
The latter does not contribute to the width at zero angle scattering in the horizontal plane 
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FIG. 2: Correlation of x-coordinates as measured by the two detectors of arm A before cuts were 
applied. Note that the background appears enhanced due to the saturation of the main band. 




FIG. 3: The raw asymmetry £{(j)) for the full |t|-interval. 

but contributes significantly at large scattering angles. Thus the x-coordinate of correlation 
distribution with the minimal width defines the position of tracks scattered at zero angle (or 
the position of the beam in the detectors) in the horizontal plane. The mean coordinates 
of non-scattered beams in the detectors were used for planar scattering angle determination 
instead of the mean coordinates of IP and mean beam angles. This approach eliminates the 
contribution of the detector position survey errors, see also discussion of systematic errors. 
To select an elastic event, a match of hit coordinates (x,y) from detectors on the opposite 
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sides of the IP was required to be within 3a for x and y-coordinate. The hit coordinates 
(x,y) of the candidate proton pairs were also required to be in the acceptance area of the 
detector, determined by the aperture of the focusing quadrupoles located between IP and 
the RP's. In case that there were more than one match between the hits on opposite sides 
of the IP the following algorithm was applied. If there is only one match with number of 
hits equal to 4, it is considered to be the elastic event. If there is no match with 4 hits or 
there are more than one such match, the event is rejected. 

The average detector efficiency was 0.98, and the upper bound of the elastic events loss 
due to all criteria was < 3.5%. 

The background originates from particles from inelastic interactions, beam halo particles 
and products of beam-gas interactions. The estimated background fraction varies from 
0.5% to 9% depending on the y-coordinate. Since in our analysis the coordinate area was 
essentially limited to ?/ > 30 strips, the background in the final sample does not exceed 2%. 



IV. DETERMINATION OF ANALYZING POWER An 

After the above cuts, the sample of 1.14 million events, for A^^^ and N^^ bunch com- 
binations, in the t-interval 0.010 < —t < 0.030, subdivided into three intervals 0.010 < 
-t < 0.015, 0.015 < -t < 0.020, 0.020 < -t < 0.030, was used to determine An- In each 
t-interval the asymmetry was calculated as a function of azimuthal angle using 5°-bins. 
Azimuthal angle dependence of the cross section for the elastic collision of the vertically 
polarized protons is given by 

^""1771 = ^-('^ + (Pb + Py)An COS0 + PbPy{Ann cos^ + Ass sin' 4>)) , (2) 
dtdcp dt 

where Pb and Py are the beam polarizations and A^^, Ass are double spin asymmetries (see 
Ref. 0] for definitions). Then the square root formula |15] for the single spin raw asymmetry 
e{(f)) can be written as 

{Pb + Py)An cos(f) 



1 + PbPy{Ann cos2 + Ass sin' 



v/A^n(0)ArU(vr-0)- ^A^U(0)Arn(7r -0) 
V^TT(0)iVU(p^ 0) + ^A^ii(0)A^TT(vr - 0) ' 
Beam polarizations for our run were |l6|Py = 0.345 ± 0.066 and Pb = 0.532 ± 0.106, 
leading to an upper constraint of 0.028 for the term Pb Py {A n n cos^ (p + Ass siia^ (p), even 



if both double-spin asymmetries Ann and Ass were as large as 0.15. This term is small 
in comparison to the systematic errors on An and was therefore neglected in Eq. Q but 
included in the systematic error, as described below. A cosine fit to the raw asymmetry e{(j)) 
was used to determine values of An, see Fig. |21 

V. SYSTEMATIC ERRORS 

Equation (jH)), from which the asymmetry is calculated has important features; namely, 
luminosities of the differently polarized proton beam bunches cancel as do the relative de- 
tection efficiencies, including geometrical acceptance, for each t and (p. 

However, two other contributions to the systematic error have to be considered: back- 
grounds, which affect the asymmetry value, and sensitivity to the transport matrix param- 
eters and to the beam position with respect to the detectors that affect the determination 
of t and 0. 

To check the effect of background, additional selection criteria were applied: 1) rejection 
of the events with a hit in one of the two y-strips closest to the beam; 2) rejection of events 
close to the boundary in the plane. From these studies, we have found that the upper 
limit of the systematic error due to the background is 4.5%. 

The final results were obtained with a transport matrix, which was obtained by correct- 
ing the standard transport matrix provided to us by the accelerator physicists from the 
Collider-Accelerator Department (C-AD). The corrections were calculated using the fully 
reconstructed tracks in all four RPs. The results were compared with those obtained with 
the standard transport matrix. The relative difference in An for the two cases is 1.4%. The 
systematic error due to an uncertainty of beam positions at the detectors is 1.8%. 

Sensitivity to the variation in Lg// was also studied and estimated to be 6.4% assuming 
upper values of transport uncertainties of Lgj^ and L^j^j as large as AL^jj-fL^j^jr = 0.1 and 
ALljrj/ Lljj- = 0.05, correspondingly. 

As mentioned earlier, neglecting the term with double-spin asymmetries in formula Q 
results in an error 2.8%. 

Since all the above errors are uncorrelated adding them in quadrature results in the 
systematic error of AAn/An = 8.4%. This error is smaller than the statistical errors of the 
measurement, cf. Table HI 
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The polarization values of the proton beams were obtained from the C-AD [l6[. They 
were evaluated using measurements for elastic proton-Carbon (pC) scattering at small 
|t|-values, in the range 0.01 — 0.02 (GeV/c)^. The details are described in Ref. |lO|. During 
the period in 2003 when the present data were taken the beam polarizations were Py = 
0.345 ±0.066 and Pb = 0.532 ±0.106. The errors include the contribution of the systematic 
part of the error due to the calibration of pC polarimeter of 13%, which is correlated for both 
beams and the uncorrelated statistical errors of the measurement. This gives the statistical 
and systematical errors of the measurement in the sum of the polarizations Py ± Pb = 
0.877 ±0.149. 

The total systematic error is comprised of A^r scale error of 17.0% mostly due to the 
systematic error of the polarization measurement, and 8.4% error due to the experimental 
systematic effects as described above. 

An important check of a possible false asymmetry e' was obtained from the asymmetry 
calculated for spin combinations A^^^ and A^^^ with a formula similar to 0. This term 
is e' ~ {Pb — Py) ■ A^. Given that the polarization values for N^^ and A^-'-^ bunches were 
Py = 0.476 ±0.085 and = 0.430 ±0.089 and the An in our t-range, one gets e' = -0.0011 
to be compared with the value we measured -0.0016, a good agreement indicating that there 
is no major source of a false asymmetry. 



TABLE 1: An results 



-t interval (GeV/c)^ 


0.010-0.015 


0.015-0.020 


0.020-0.030 


0.010-0.030 


<-t> (GeV/c)2 


0.0127 


0.0175 


0.0236 


0.0185 


An 


0.0277 


0.0250 


0.0178 


0.0212 


A An - stat. 


±0.0061 


±0.0043 


±0.0030 


±0.0023 


AAn - syst.* 


±0.0023 


±0.0021 


±0.0015 


±0.0018 


AAn due to A{Py + Pb) 


±17.0 % 


* Contributions to systematic error were added in quadrature 
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FIG. 4: The single spin analyzing power for three t intervals. Vertical error bars show statistical 
errors. The solid curve corresponds to theoretical calculations without hadronic spin-flip and the 
dashed one represents the fit. 

VI. RESULTS AND CONCLUSIONS 

The values of obtained in this experiment and their statistical errors are shown in 
Fig. 0] for the three t-intervals, and they are summarized in Table IH 

The curves shown in the figure represent theoretical calculations using the formula for 
An in the CNI region. The general formula is given by Eq. 28 of Ref. 0]. With reasonable 
assumptions that the amplitude 02 and the difference 0i — 03 could be neglected at collider 
energies, the formula becomes simpler 

4 [k{1 - p5) + 2{6 Re - Im rs)]^ - 2(Re rs - p Im r^) 



A 



N 



(4) 



m (k)2_2(p + 5)k + (i + p2) 

In this formula tc = —Sna/atot, is the anomalous magnetic moment of the proton, p is 
the ratio of the real to imaginary parts of forward (nonflip) elastic amplitude, and S is the 
relative phase between the Coulomb and hadronic amplitudes. Since the total cross section 
(o'tot) and the p parameter have not been measured in this energy range, we have used values 
of atot = 51.6 mb and p = 0.13. These values come from fits to the existing pp data taken 
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FIG. 5: Fitted values of (full circle) with contours corresponding to the different confidence 
levels. The point corresponding to no hadronic spin-flip (triangle) is also shown. 

at energies below 63 GeV and world pp data. They also agree well with the predictions of 



various models 



Id 



The Coulomb phase 6 is calculated as in Ref. 
2 



5 = a In 



a 7, 



(5) 



|t|(6 + 8/A2) 

where b is the slope of the forward peak in elastic scattering, a is the fine structure constant, 
Euler's constant 7 = 0.5772 and = 0.71 GeV^. The value of b comes from our previous 
measurement 

The solid curve in Fig. |3] corresponds to the calculation without hadronic spin- flip (Re rs 
and Im set to in Eq. 0]). To quantify a possible contribution of the single helicity- 
flip amplitude 05, the formula given by Eq. I^was fitted to the measured Aj^ values with 
Re rs and Im as fit parameters. The statistical and systematical errors (except the beam 
polarization error) of An were added in quadrature for the fit. The results of the fit are 
following: Re = —0.033 ± 0.035 and Im = —0.43 ± 0.56. The dashed line in Fig. |3] 
respresents the curve resulting from the fit. 

The fitted values of Rers and Imrs are shown in Fig. Altogether with contours for Icr, 2a 
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and 3(T confidence levels. In addition, the point corresponding to no hadronic spin-flip is also 
shown. The fitted rs is compatible, at about one a level, with the hypothesis of no hadronic 
spin flip. Thus our conclusion is that our results are suggestive of a hadronic spin-flip term, 
but cannot deflnitively rule out the hypothesis that only hadronic non spin-flip amplitudes 
contribute. 

Recent measurements of Ajv at substantially lower cms energies than the one reported 
here indicate small, but signiflcantly different from zero, contribution of spin-flip amplitude 
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in case of proton-carbon scattering 
for proton-proton scattering ^ at ^/s 



and are consistent with no spin-flip contribution 
13.7 GeV. 
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